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ABSTRACT 


Single-layer mixed-phase stratiform (MPS) Arctic clouds, which formed under conditions 
of large surface heat flux combined with general subsidence during a subperiod of the 
Atmospheric Radiation Measurement (ARM) Program Mixed-Phase Arctic Cloud Experiment 
(M-PACE), are simulated with a cloud resolving model (CRM). The CRM is implemented with 
either an advanced two-moment (M05) or a commonly used one-moment (L83) bulk 
microphysics scheme and a state-of-the-art radiative transfer scheme. 

The CONTROL simulation, that uses the M05 scheme and observed aerosol size 
distribution and ice nulei (IN) number concentration, reproduces the magnitudes and vertical 
structures of cloud liquid water content (LWC), total ice water content (IWC), number 
concentration and effective radius of cloud droplets as suggested by the M-PACE observations. It 
underestimates ice crystal number concentrations by an order of magnitude and overestimates 
effective radius of ice crystals by a factor of 2-3. The OneM experiment, that uses the L83 
scheme, produces values of liquid water path (LWP) and ice plus snow water path (ISWP) that 
were about 30% and 4 times, respectively, of those produced by the CONTROL. Its vertical 
profile of IWC exhibits a bimodal distribution in contrast to the constant distribution of IWC 
produced in the CONTROL and observations. 

A sensitivity test that uses the same ice-water saturation adjustment scheme as in OneM 
produces cloud properties that are more similar to the OneM than the CONTROL. The 
CONTROL predicts spatially varying values of the intercept parameter of snow size spectra (N 0s ) 
that are one order of magnitude smaller than the prescribed N 0s used in L83. A sensitivity test that 
prescribes the larger L83 N 0s results in 20% less LWP and 5 times larger snow water path than the 
CONTROL. When an exponential ice size distribution replaces the gamma size distribution in the 
CONTROL, ISWP decreases by 70% but LWP increases by 7% versus the CONTROL. 

Increasing the IN number concentration from the observed value of 0.16 Z' 1 to 3.2 L~ ] forces the 
MPS clouds to become glaciated and dissipate, but the simulated ice number concentration agrees 
initially with the observations better. Physical explanations for these quantitative differences are 



provided. It is further shown that the differences between the OneM and the CONTROL are larger 
than those due to the estimated uncertainties in the prescribed surface fluxes. Additional 
observations and simulations of a variety of cases is required to further narrow down uncertainties 
in the microphysics schemes. 



1. Introduction 


Atmospheric numerical models with a horizontal grid spacing of 1 - 2 km are known as 
cloud-resolving models (CRMs). CRMs are able to resolve convective-scale and mesoscale 
circulations and, hence, can better represent the interactions between physical processes involving 
smaller scales than traditional General Circulation Models (GCMs). Physical processes such as 
those involving clouds and precipitation cannot be explicitly resolved and have to be 
parameterized in GCMs because of the grid spacing of a GCM, which is typically on the order of 
100 km in the horizontal and 1 km in the vertical. Unfortunately, there are large uncertainties in 
parameterizations of subgrid scale processes and improvement of parameterizations has been 
slow in spite of the enormous efforts made over the past decades (Randall et al. 2003). Moreover, 
the subgrid-scale processes interact mainly through the time-varying large-scale variables (and 
surface conditions) in GCMs while in reality they directly interact with each other. A unified 
formulation of the entire spectrum of these interactions is necessary for more accurate climate and 
weather prediction, but it is difficult to achieve this with the traditional grid spacings used in 
GCMs (Arakawa 2004). Therefore, with the rapid growth of computational capacity, continental- 
scale NWP is currently performed at cloud-resolving scales (e.g. WRF; Skamarock et al. 2005). 
For climate simulation, CRMs have been used as a “super-parameterization” to replace most of 
the traditional parameterizations in each grid cell of GCMs (e.g. Grabowski 2003) and global 
versions of CRMs are emerging (Tomita et al. 2005). 

Microphysical processes, as well as turbulent and radiative transfer processes, still need to 
be parameterized in CRMs. Most CRMs rely on bulk microphysics schemes to represent the 
complicated interactions between atmospheric thermodynamic states and hydrometeors and 
among various hydrometeor species. Bulk microphysics schemes typically divide the 
hydrometeor spectrum into cloud water, cloud ice, rain, and one or more ice-phase precipitation 
species (e.g. snow, graupel, and hail). Each hydrometeor class is represented by a specified size 
distribution function (e.g. gamma, exponential, and lognormal). The microphysics schemes that 
predict only hydrometeor mixing ratios are called the one-moment approach (e.g. Lin et al. 1983). 
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An improvement to the one-moment approach is to predict the rates of change for both mixing 
ratios and number concentrations of hydrometeors, i.e. the two-moment approach (e.g. Ferrier 
1994; Meyers et al. 1997; Morrison et al. 2005, hereafter M05; Vaughan et al. 2007). An 
advantage of this approach is that the effective sizes of cloud particles, one of the most important 
parameters determining cloud radiative impacts, can be predicted, in contrast to the one-moment 
approach. Another advantage is that two-moment schemes potentially can represent the size 
distributions of hydrometeors more realistically and thus represent microphysical processes more 
accurately than one-moment schemes (e.g., Meyers et al. 1997; Morrison and Pinto 2006). 

Arctic clouds have been identified as playing a central role in the Arctic climate system. 
However, the role of clouds is even less well understood in the Arctic than in other geographic 
regions, due to sparse observations. The Arctic field programs such as the Surface Heat Budget of 
the Arctic (SHEBA; Uttal et al. 2002) and the First ISCCP Regional Experiment (FIRE; where 
ISCCP is the International Satellite Cloud Climatology Program) Arctic Cloud Experiment (ACE; 
Curry et al. 2000) revealed that mixed-phase stratiform (MPS) clouds appear to dominate the low- 
cloud population within the Arctic (Intrieri et al. 2002). Moreover, it is found that the Arctic 
mixed-phase clouds are distinct from their lower latitude cousins (e.g. Curry et al. 1996, 2000). A 
unique feature of these clouds is that they are persistent, liquid-topped clouds that precipitate ice 
(Hobbs and Rangno 1998; Intrieri et al. 2002). Another unique feature of these clouds is that the 
liquid component of the mixed-phase cloud dominates the radiative properties (McFarquhar and 
Cober 2004; Zuidema et al. 2005) 

Adequate simulation of Arctic clouds is needed to address Arctic cloud-radiative-surface 
interactions that may impact global climate (e.g. Curry et al. 1996) and to predict weather, due to 
the persistence and large horizontal extent of these cloud systems. However, there have been few 
simulations of Arctic MPS clouds with CRMs, primarily because the observations of cloud 
physical properties needed to evaluate model performance are sparse and there is a lack of large- 
scale forcing data available to drive CRMs. The Department of Energy - Atmospheric Radiation 
Measurement (DOE- ARM) Program (Stokes and Schwartz 1994; Ackerman and Stokes 2003) 
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recently launched its Mixed-Phase Arctic Cloud Experiment (M-PACE; September 27 - October 
22, 2004) at the North Slope of Alaska (NS A) sites (Harrington and Verlinde 2004; Verlinde et al. 
2007). During the field campaign, detailed information about Arctic clouds were measured using 
the ARM millimeter wave cloud radar, micropulse lidar, laser ceilometers, microwave radiometer 
(MWR), and three instrumented aircraft. Furthermore, the large-scale forcing data were derived 
for a seventeen and a half day Intensive Observation Period in October 2004 (Xie et al. 2006) by 
applying the method of Zhang and Lin (1997) and Zhang et al. (2001) to the available data. These 
forcing data can be used to drive models [CRMs, single-columns models (SCM; Randall et al. 
1996), and large-eddy simulation (LES) models]. 

The objectives of this study are two-fold. One is to evaluate CRM simulations of Arctic 
MPS clouds with a state-of-the-art dataset. The available M-PACE data offer a promising 
opportunity for improving cloud microphysical parameterizations in CRMs. Here, single-layer 
MPS clouds observed during a sub-period of M-PACE are simulated using a CRM, driven by the 
ARM-derived large-scale forcing. The CRM includes a state-of-the-art radiative transfer scheme 
and either a one- or a two-moment microphysics scheme. The performance of the CRM is 
evaluated through comparing simulated cloud properties, such as the vertical profiles of cloud 
liquid water content (LWC), ice water content (IWC), droplet number concentration, ice number 
concentration, effective sizes of droplets and ice crystals, with the M-PACE aircraft observations 
(McFarquhar et al. 2007), as well as the retrievals of liquid water path (LWP; Turner et al. 2007) 
and observations of precipitation from ground-based instruments deployed at the NS A sites. 

The second objective of this study is to explore the sensitivities of the simulated clouds to 
representation of various microphysical processes and parameters. To achieve this objective, a 
range of sensitivity tests are conducted. In particular, we attempt to answer the following 
questions: what differences in the simulated cloud properties are produced by use of a one- or 
two-moment microphysics approach? What microphysical processes and parameters may 
significantly influence the simulated MPS clouds? 
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Section 2 describes the CRM used in this study, with a focus on the prediction of 
hydrometeor number concentrations. Section 3 gives a description of the case and cloud-property 
observations. Design of the numerical experiments is presented in Section 4. Results from the 
CRM simulations utilizing either the one- or two-moment approach are compared with the 
aircraft measurements in Section 5. Section 6 contains results from the sensitivity tests. Summary 
and conclusions are given in Section 7. 

2. The numerical model 

The dynamic framework of the CRM used in this study is based on the anelastic forms of 
hydrostatic, momentum and continuity equations in two dimensions (x andz) with a third-moment 
turbulence closure (Krueger 1988; Xu and Krueger 1991). The CRM includes the Fu-Liou (1993) 
radiative transfer parameterization and either a one-moment or a two-moment microphysics 
parameterization. The two-moment bulk microphysics scheme of M05 has been implemented, 
which predicts the mixing ratios and number concentrations of cloud water, cloud ice, rain, and 
snow. The equations used to predict the hydrometeor number concentrations are: 


dn^ 

dt 


-^! ( Po <»")+ a , + s x + m x 


( 1 ) 


where n x is the number concentration with the subscript x being c, i, r, s for cloud water, cloud 
ice, rain, and snow, respectively. p 0 is the dry air density of the initial (reference) state, n” w" is 
the ensemble mean of the turbulent flux of n x in the vertical direction. A x refers to activation (for 
cloud water) and nucleation (for cloud ice), S x represents sedimentation, and M x denotes all other 
microphysical processes. The effects of turbulent fluctuations on number concentrations of 
raindrops and snow are ignored in the current version of the CRM, and the effects of turbulence 
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on their mixing ratios are also ignored (Krueger 1988). For number concentrations of cloud water 
and ice, K-theory is applied to determine the turbulence terms; that is, 

Po<^ = -Po*^T ( 2 ) 

The exchange coefficient K is calculated using K = cl J TKE , where c is a constant (0.24), / is 
the turbulence length scale and TKE is the turbulent kinetic energy. Both / and TKE are 
determined by the third-moment turbulence closure (Krueger 1988). 

Droplet activation is treated by a physically -based scheme (Abdul-Razzak et al. 1998; 
Abdul-Razzak and Ghan 2000). This scheme not only relates droplet activation to aerosol 
characteristics but also couples it with local cooling rate that is determined by cloud-scale and 
sub-grid turbulent vertical velocity as well as radiative cooling. The error of the parameterization 
is less than 10% under a wide variety of conditions (Abdul-Razzak et al. 1998; Abdul-Razzak and 
Ghan 2000). The turbulent upward motion for the droplet activation calculation is approximated 
as the square root of the vertical component of TKE per unit mass. Sedimentation of cloud 
particles is calculated with terminal particle fall velocities related to particle sizes and air density 
(Ikawa and Saito 1991). Parameterizations of all other microphysical processes follow M05, 
including deposition, condensation-freezing of ice nuclei, contact- and immersion-freezing 
nucleation of cloud droplet and raindrops, autoconversion of cloud water to rain and of cloud ice 
to snow, self-collection of cloud droplets and of raindrops, snow aggregation, accretion of cloud 
droplets, rain and cloud ice by snow, rime-splintering from accreted droplets and raindrops by 
snow, accretion of cloud water by rain, deposition/sublimation of cloud ice and snow, melting of 
snow, evaporation of rain and melted snow, saturation adjustment of cloud water, as well as the 
decrease in number concentrations during evaporation/sublimation. 
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In the M05 scheme, the gamma size distribution is assumed for cloud droplets and cloud 
ice crystals while the Marshall-Palmer (exponential) size distribution is used for raindrops and 
snowflakes. A gamma size distribution can be expressed as 


N(D) = N 0 D^e XD 


(3) 


where D is diameter, N 0 is the “intercept” parameter, p is the spectral shape parameter, and X is 
the slope parameter. The value p is determined by the relative radius dispersion (p; defined as the 
ratio between the standard deviation and the mean radius): 

n = i /n 2 - 1 (4) 

Practically, parameters N 0 and X can be diagnosed from the specified p and predicted mixing ratio 
( q ) and number concentration (n) of the species. That is, only p needs to be specified using the 
two-moment approach. For the one-moment approach, two of the three parameters (N 0 , p, and X) 
need to be specified. 

For cloud droplets, p is related to the number concentration, n c , in the M05 scheme. 
However, the exact p -n c relationships for Arctic clouds are not yet developed. There are currently 
only a few formulations relating p to n c and these are based on observations at lower latitudes. For 
example, Rotstayn and Liu (2003; RL03) fitted three curves to measurements in polluted and 
unpolluted warm stratiform and shallow cumulus clouds. These curves are designed to represent 
the average variation of p with n c , as well as lower and upper bounds of this variation. These 
curves shown in Fig. la are defined by 
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r) = 1 -0.7 e 


( 5 ) 


-an r 


where a equals 0.001 for the lower curve, 0.003 for the middle curve, and 0.008 for the upper 
curve. The corresponding \i-n c relationships are displayed in Fig. lb. The relationship of Eq. (5) 
with a of 0.003 is used in this study. 1 * Note that there was considerable scatter in the data used by 
Rotstayn and Liu (2003) to obtain the r\-n c relationship of Eq. (5). Miles et al. (2000) created a 
database of stratus cloud droplet size distribution parameters, derived from in-situ data reported in 
the existing literature. The datasets included several parameters for 42 marine stratocumulus 
clouds and 52 continental stratocumulus clouds. These observations, however, do not show a 
systematic increase or decrease in r\ with increasing n c . For cloud ice, a constant p of 5 is used in 
M05, corresponding to a q of -0.408. Note that the Marshall-Palmer distribution is a special case 
of Eq. (4) with p equal to zero. For the radiation calculation, the effective sizes of cloud water, 
cloud ice and snow are determined by the predicted size distributions. 

The CRM also includes the commonly used one-moment bulk microphysics scheme of 
Lin et al. (1983) (L83 hereafter) with modifications to its ice-phase microphysics 
parameterization by Krueger et al. (1995). This scheme represents the rates of change of mixing 
ratios for five hydrometeor species (cloud water, cloud ice, rain, snow, and graupel) and is 
combined with an ice-water saturation adjustment (Lord et al. 1984) to determine the 
condensation/evaporation of cloud water and deposition/sublimation of cloud ice. Cloud water 
and cloud ice are assumed to be monodisperse. Precipitating hydrometeor species are assumed to 
have exponential size spectra. Number concentrations of the precipitating hydrometeor species 
can be diagnosed from the predicted mixing ratios and specified microphysical parameters 

1 We also tested formulations for the spectral shape parameter (p) as a function of n c that were used in Grabowski 

(1998) and Morrison and Grabowski (2007). Figure 1 indicates that these formulations produce substantially dif- 

ferent r| at most values of n c . For values of n c of ~60 cm" 3 simulated for this case study, however, the results are 
not sensitive to the specific formulation of p. Therefore, these results are not included in this paper. 
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describing the hydrometeor size spectra. However, aerosol characterization is not physically 
linked to the hydrometeor number concentrations. For the radiation calculation, the effective 
radius of cloud droplet is specified (10 pm) and the effective sizes of cloud ice and snow are either 
empirically determined from IWC or specified (120 pm for snow), as in our earlier studies (Xu 
2005; Luo et al. 2007). 

The water-ice saturation adjustment scheme of Lord et al. (1984) requires assumptions 

about both the coexistence of cloud water and cloud ice at temperatures less than 0°C and the 
partitioning between condensation and deposition. Specifically, the Lord et al. scheme assumes 
that the saturation vapor mixing ratio q* is a mass-weighted average of the respective saturation 

values over liquid water and ice at -40 °C <T<0°C when both cloud water and cloud ice are 
present. Under subsaturated conditions, cloud water is evaporated first so that water vapor mixing 
ratio (q v ) would be equal to q*. If subsaturated conditions are still present after all cloud water 

evaporates, enough cloud ice is sublimated such that q v <q*. On the other hand, production of 

either cloud water (Aq c ) or cloud ice (A q t ) depends linearly on temperature under supersaturated 

conditions so that A q c = q v - q* at T = 0°C and A q t = q v -q* at T = -40 °C . A similar 

formulation was also developed by Tao et al. (1989) except for removing the iterative adjustment 
procedure used in Lord et al. (1984). 

3. Description of the case study 

The east-northeast flow brought cold near-surface air from the sea-ice located about 500 
km north over the warm open ocean that was adjacent to the northern coast of Alaska (Fig. 2). The 
contrast between the cold-air and warm open ocean resulted in large ocean sensible and latent 
heat fluxes which, combined with the conditions of large-scale subsidence, promoted a well- 
mixed cloudy boundary layer. Single layer mixed-phase clouds were formed under these 
conditions (Verlinde et al. 2007). These clouds were then advected to the Alaskan coast where 
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they were observed at the ARM NSA sites — Barrow and Oliktok Point (Fig. 2). The ARSCL 
(Active Remote Sensing of Clouds) algorithm (Clothiaux et al. 2000) derived cloud distribution 
exhibits the presence of single layer stratocumulus in the period 9-14 October, 2004 (not shown). 
The time-height distributions of radar reflectivity, lidar backscatter and lidar depolarization (e.g. 
Figure 6 of Verlinde et al. 2007) reveal the locations of cloud top and cloud liquid base, and the 
presence of shafts of ice precipitation and/or drizzle throughout the cloud layer and below cloud. 

The bulk microphysical properties of the MPS clouds that occurred during M-PACE, i.e. 
total condensed water content, LWC, IWC, effective radius of supercooled water droplets, 
effective radius of ice crystals [defined following Fu (1996)], total water droplet number 
concentration and total ice crystal number concentration, were derived by McFarquhar et al. 
(2007) from measurements obtained by instruments on the University of North Dakota Citation 
aircraft. The Citation was equipped with a range of probes for measuring the size, shapes, and 
phases of the complete range of hydrometeors that can be sampled within a cloud. There were one 
Citation flights on October 9 and 12, respectively, and two on October 10, which occurred in 
single layer MPS clouds that were similar in structure. The four flights covered a period of ~ 6.5 h 
with about half of the period for in-cloud observation. Here the cloud base is defined as the lidar- 
derived liquid cloud bottom. The cloud top is defined as the cloud radar-derived cloud top or, 
when cloud radar data was not available, as the location where the total condensed water content 

became greater than 0.001 g m' 3 (McFarquhar et al. 2007). The bulk properties are available at 10 
s resolution, but represent a 30 s running average of the measured ice properties. There are 1131 
in-cloud samples obtained from the four flights. The bulk cloud properties sampled by the four 
flights are used to validate model simulations in this study. 

Other evaluation data include measurements of LWP provided by the microwave 
radiometer (MWR) (Turner et al. 2007) and those of surface precipitation provided at Barrow site. 
Large uncertainties, however, existed in the ARM surface precipitation measurements during M- 
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PACE because of both the blowing snow conditions and the lack of a dense observational network 
(Xie et al. 2006). 

4. Design of CRM simulations 

We conduct a set of simulations using the CRM described in Section 2 to explore the 
model ability to simulate the MPS clouds and its sensitivity to microphysics scheme and 
parameter. All these simulations start with the same initial profiles of the atmospheric state. They 
are prescribed with the same surface latent and sensible fluxes, large-scale subsidence, and 
horizontal advection of temperature and moisture. Details of forcing data are described in Section 
4a. For the sensitivity simulations, different treatments of some microphysical processes and 
parameters, described in Section 4b, are used. The horizontal grid spacing is 2 km. The vertical 
grid spacing varies with height from 30 m to 102 m at heights below 1 .9 km and is constant (500 
m) above 1 .9 km. The domain width is 256 km in the horizontal and 20 km in the vertical. A time 
step of 5 seconds is used for all simulations. 

4.1 Initial conditions, large-scale forcing, and aerosol specification 

The initial and lower boundary conditions, large-scale forcing data, and aerosol properties 
provided by Klein et al. (2006) are used in all simulations. The period of our simulation is from 
17Z October 9 to 5Z October 10. The initial profiles of temperature and water vapor are based on 
the 17Z October 9 sounding at Barrow (Figs. 3a, b) with the inversion height at ~1.4 km. The 
CRM is initialized with an adiabatic profile of liquid water (Fig. 3b). No ice is present at the initial 

time. The total water mixing ratio below inversion is 1.95 g kg' 1 . The CRM starts from 
horizontally homogeneous fields except for the added random perturbations with a maximum of 
0.1 K to the potential temperature field at the lowest several levels. 

The forcing data were based on an analysis of the ECMWF model data for the oceanic 
region adjacent to the NSA sites (Xie et al. 2006). The magnitude of the large-scale subsidence 
(&>) linearly increases with decreasing pressure from a zero value at the surface to a value of about 
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3.3 hPa h" 1 at and above the inversion (Fig. 3c). This is used to vertically advect all 
thermodynamic and microphysical variables in the model. The large-scale horizontal advective 
tendencies of temperature and moisture are also prescribed (Klein et al. 2006; also shown in Figs. 
3d, e). Due to the lack of observations, the large-scale horizontal advective tendency of the cloud 
variables are set to zero. The CRM’s horizontally-averaged winds (u and v) are also nudged 

toward the initial values (-13 m s' 1 for u and -3 ms" 1 for v, respectively) with a time scale of 1 h 
(Xu and Randall 1996). Surface sensible and latent heat fluxes are specified as 136.5 W m" 2 and 

107.7 W m" 2 , respectively. For radiation purposes, the lower boundary is an open-ocean surface. 
An SST of 274.01 K is used in the upward longwave radiation calculation. The spectral surface 
albedos for the six bands of Fu and Liou (1993) radiation code are calculated using the 
parameterization of Jin et al. (2004). 

The CRM’s droplet activation parameterization is physically linked to the characterization 
of aerosols. We use a bimodal lognormal size distribution of dry aerosol, obtained lfom a Met One 
Fland-Held Particle Counter (HHPC-6) on board the ARM unmanned aerial vehicle (UAV) and a 
condensation nuclei counter from the NOAA Earth System Research Laboratory located near 
Barrow, AK. The size distribution for each mode of the lognormal distribution is represented by 


dN N t f ln 2 ( r/r m)) 

^ lnr Torino l 21n 2 <j ) 


( 6 ) 


where the parameters N t , a , and r m are the total number concentration, standard deviation, and 
geometric mean radius of each mode, respectively. For the smaller mode, the values of these 
variables are 72.2 cm' 3 , 2.04, and 0.052 pm, respectively. The corresponding values for the larger 
mode are 1.8 cm' 3 , 2.5, and 1.3 pm. The aerosol composition is assumed to be ammonium 
bisulfate with an insoluble fraction of 30%, as recommended by Klein et al. (2006) based on 
observations (Bigg and Leek 2001; Zhou et al. 2001). 
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In-situ out-of-cloud observations for number concentration of active ice forming nuclei 
(IFN) were obtained on October 9 and 10 from the Continuous Flow Diffusion Chamber (Rogers 
et al. 2001) aboard the Citation aircraft. These measurements represent the total number 
concentration of active IFN that have diameters less than 2 pm acting in deposition, condensation- 
freezing, and immersion-freezing modes. The measured mean concentration of these IFN is about 

0.16 L' 1 , which is used to represent the aforementioned nucleation modes in the CRM 
simulations. 

4.2 Sensitivity tests 

In order to explore the possible impacts of microphysical processes and parameters on 
CRM-simulated MPS clouds, a range of sensitivity tests are performed (Table 1). The baseline 
simulation (hereafter referred to as CONTROL) is performed with a two-moment approach for 
both cloud particles and precipitating hydrometeor species using the M05 scheme. A sensitivity 
experiment, OneM, is performed with a one-moment approach for all hydrometeor species as 
described in Section 2 to quantify the benefits of the two-moment approach. Note that graupel is 
allowed to occur in the OneM simulation but it never does. 

A sensitivity test (SAT), which is the same as the CONTROL except for using the water- 
ice saturation adjustment scheme of Lord et al. (1984), is designed to examine the role of the 
water-ice saturation adjustment used in the one-moment microphysics parameterization (Lord et 
al. 1984; Tao et al. 1989). The Lord adjustment scheme, described in Section 2, is different than 
the M05 scheme used in the CONTROL, which determines deposition/sublimation of cloud ice 
(as well as snow and rain) using a non-steady, vapor diffusion approach and applies a saturation 
adjustment approach only to cloud liquid water, which is reasonable because of short droplet 
phase relaxation time. 

The rest of microphysics experiments test several microphysical parameters used in the 
M05 scheme. Experiment IN20 is performed by increasing the IFN number concentration by a 

factor of 20 from the measured value, i.e. from 0.16 L" 1 to 3.2 L" 1 . This experiment is motivated 
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by previous numerical modeling studies of Arctic MPS clouds which showed large sensitivity of 
simulated MPS clouds to the availability of IFN (Harrington et al. 1999; Jiang et al. 2000; 
Morrison and Pinto 2006). Morrison and Pinto (2006) found that the prediction of n s could 
critically affect an MPS cloud simulated by a mesoscale numerical model. To examine this issue, 
experiment NOS is performed by setting the intercept parameter N 0s equal to a constant value of 

3.0E6 m' 4 (Gunn and Marshall 1958; Lin et al. 1983) so that the number concentration of snow 
particles, n s , is diagnosed rather than predicted. The last sensitivity test, piO, examines the spectral 
shape parameter (p) in the gamma size distribution (Eq. 4) of cloud ice in the two-moment 
approach. Experiment piO is performed with pj of zero, instead of 5 in the CONTORL. That is, 
cloud ice is represented by an exponential (rather than a gamma) size distribution in the piO 
experiment. 

Another set of sensitivity tests (Table 1) aim at examining the impacts of estimated 
uncertainties in the surface fluxes, which are compared to the differences between the one- 
moment and two-moment schemes. These tests are the same as either the CONTROL or the 
OneM simulations, except for increasing or decreasing the surface sensible and latent heat fluxes, 
respectively, by 10%. One reason for performing these tests is that the magnitudes of these fluxes 
were based on the ECMWF model data for the oceanic region adjacent to the NSA sites and, 
therefore, may contain model uncertainties. Another reason is that previous modeling studies 
indicate that surface turbulent flux could influence properties of simulated mixed-phase Arctic 
clouds (e.g., Harrington and Olsson 2001). 

5. Comparison between CRM simulations and aircraft observations 

We first examine the CONTROL and OneM simulations since they represent results using 
the two distinct (two-moment vs. one-moment) microphysics schemes. 

5.1 Vertical profiles of hydrometeor mass 
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The vertical profiles of LWC and IWC plus snow water content (hereafter, ISWC) from 
the CONTROL and OneM simulations and observations (means plus/minus standard deviations 
computed from the four flights) are compared to examine the vertical variations of cloud 
distributions. The simulated LWC and ISWC are horizontally averaged and time averaged at 30 
min blocks during the 12 h simulation period. Only those centered at 3.25 h, 10.25 h, and 1 1 .75 h 
are shown in Fig. 4. The observations represent both spatial and temporal variability since many 
of the observations were obtained in different locations (Barrow, Oliktok Point and in between). 
Following McFarquhar et al. (2007), the vertical axis of Fig. 4 is a normalized height (H n ) defined 

as (H - H b ) /(H t - H b ) , where H is the height, H b cloud base height and H t cloud top height. 

The cloud top and cloud base are located at H n = 1 and H n = 0 , respectively. A negative H n 

represents a height below the liquid cloud base. Observations below liquid cloud base typically 
refer to the presence of precipitating ice, and on occasion refer to an erroneously identified cloud 
base. The observations are categorized into 20 bins of H n within the cloud layer. There are about 
50 samples for each of the observed cloud properties within each H n bin. 

McFarquhar et al. (2007) analyzed the variation of the observed microphysical variables 
with height. In order to compare against the model simulations, the most notable features are 
summarized here. The observed, averaged LWCs increase with height within the cloud layer with 
a peak of ~0.32 g m' 3 located near the cloud top. The standard deviations of the observed LWC 
range from 0.05 g m" 3 to 0.08 g m" 3 below cloud top (H n < 0.8) and increase to ~0. 14 g m' 3 at the 
cloud top. The larger variation of the observed LWC near cloud top may be related to entrainment. 
The observations also indicate that there is a small amount of ISWC (0.01 g m" 3 ) with a relatively 

constant vertical distribution within the cloud layer, but with large variations (up to 0.04 g mf 3 ) in 
the lower part of the cloud layer (H n < 0.25). The large variations suggest that large ISWCs were 
only occasionally observed near cloud base. The observed fraction of ice to the total condensed 
water, however, increases towards the base of the cloud (McFarquhar et al. 2007). 
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The liquid and ice coexist throughout the entire period of the two simulations (Figs. 4a, 
b), consistent with the observations which showed mixed-phase clouds occurred 71% of the time 
for the observations. The cloud top and cloud base in the model are located at 1.33 km and 0.65 
km, respectively. In both simulations, ice crystals (including snow) occur throughout the cloud 
layer and fall below liquid cloud base to the surface, consistent with radar and lidar measurements 
shown in Verlinde et al (2007). However, there are some obvious differences between the two 
simulations. In the CONTROL, both the LWC and the ISWC reach a steady state after ~3 h. The 
LWCs increase with height and the ISWCs are constant with height within the cloud layer. Both 
the LWC and ISWC are located within the uncertainty range of the observations. In the OneM 
experiment, the liquid cloud layer decays with time and the ice mass increases with time. The 
amount of LWC is underestimated compared to the observations. The ISWCs from the OneM 
experiment exhibit larger variations with height as well as larger amounts at most heights within 
the cloud layer than those in the observations or in the CONTROL results. 

To further explore the differences in ice crystal mass between the CONTROL and OneM 
simulations, separate vertical profiles of IWC and snow water content (SWC) from the two 
simulations are compared (Figs. 4c, d). The IWCs from the CONTROL are nearly constant with 
height within the cloud layer. The IWCs from the OneM run exhibit two peaks, one located near 
the cloud top and the other at the lower part of the cloud layer during the majority of the 12 h 
simulation period. The only exception occurs at the last hour when there is a single peak at H n of 
~0.8. These differences are related to the cloud ice deposition process in the CONTROL and 
OneM simulations, as shown in the time-height distributions of cloud ice deposition rate in Figs. 
5a and b. In the CONTROL, deposition (from water vapor to cloud ice) occurs smoothly in height 

and in time within the cloud layer at the instantaneous rates of less than 0.01 g kg' 1 h' 1 . In the 
OneM experiment, deposition (positive values) or sublimation (negative values) rates exhibit 
significant variability within the cloud layer, and are one order of magnitude larger than those 
seen in the CONTROL. The OneM simulation also produces deposition rates (in the lower part of 
the cloud layer) that oscillate with a period of about 30 minutes during the first 8 h of the 
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simulation. After cloud water decreases significantly in the OneM simulation (i.e. after 10 h of the 
simulation), the cloud ice deposition process is enhanced significantly within the cloud layer 
because the saturation vapor mixing ratio depends on the relative amounts of liquid and ice. 

To examine the effects of microphysics on dynamics, the domain-averaged sub-grid 
turbulent kinetic energy (TKE) is compared between the CONTROL and OneM simulations 
(Figs. 5c, d). While the TKE produced by the CONTROL is relatively constant below the cloud 
top except for near the surface and changes little with time, the TKE produced in the OneM 
simulation shows more significant variability with time and height. The OneM TKE in the interior 
of the cloud layer oscillates with the same period as its deposition process does. The OneM TKE 
near the cloud base has smaller values than those in the interior of the cloud layer and below the 
cloud base. These results suggest that different representations of microphysics processes have 
distinct impacts on simulated dynamics. 

Snow exists from the cloud top to the surface in both simulations with maxima located 
near the cloud base (Fig. 4d). However, the OneM experiment produced SWCs that are several 
times as large as those from the CONTROL. The larger SWCs in the OneM experiment are partly 
attributed to the greater deposition of cloud ice (Fig. 5b) that is subsequently converted to snow 
through the autoconversion process. As will be shown in Section 6, this result can also be 
attributed to the intercept parameter of the snow size distribution (N 0s ) specified in the one- 
moment scheme, which is larger than that predicted in the CONTROL. 

5.2 Vertical profiles of hydrometeor number concentration and effective radius 

Number concentrations and effective radii of cloud liquid droplets and ice crystals are 
important cloud properties that significantly influence cloud optical properties and various 
microphysical processes. These variables are not predicted in the OneM experiment. Therefore, 
we compare those from the CONTROL simulation to the observations (Fig. 6). Averages of the 
observed droplet number concentrations (n c ) are relatively constant with height in the cloud layer 

with values of about 50 cm" 5 (Fig. 6a). The variation of the observed n c at each height bin ranges 
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from 20 to 35 cm" 3 . The simulated n c is constant with height and has a value of ~60 cm' 3 , 
generally consistent with and within the variability of the observations. The observations suggest 
that the effective radii of cloud droplets (r e ) generally increase with height within the cloud layer 
(Fig. 6b). The standard deviation of r e ranges between 1 pm and 2 pm for most height bins except 
for near the cloud top where it increased to ~ 3 pm. The simulation reproduced the observed 
increase of r e with height within the cloud layer. At most height bins, the simulated r e is within 
the uncertainty range of the observations except for near the cloud base where the observations 
are greater than the simulated r e . The underestimation of r e near cloud base is probably related to 
the smaller simulated LWC at that height, compared to the observations (Fig. 4a). On the other 
hand, on some of the ramped ascents and descents there may have been some uncertainties in the 
identification of cloud base from the measurements. 

The vertical profile of total ice crystal number concentration (n t ) from the observations 

(Fig. 6c) shows a relatively constant distribution with height with a mean of 1-3 If 1 , significantly 
greater than the observed IFN number concentration. The standard deviation of n t is comparable 
to or greater than the mean value. The simulated n i (including both cloud ice and snow) is less 

than 0.5 Z , smaller than the observed mean. Note that the observed «,• refers to concentration of 
ice particles with diameter greater than 50 pm. The discrepancy between the observed and 
simulated n t would be even larger if ice particles with diameters smaller than 50 pm were 
excluded from the simulated results. 

The observations show that the vertical profile of effective radius of ice crystals (r ei ), 
defined following Fu (1996), is constant with height and the mean values of r ei are ~ 25 pm. In the 
CONTROL simulation, ice effective radius is calculated by mass-weighting of the inverse values 
for cloud ice and snow which also follow the definition of Fu (1996). The vertical profile of 
simulated r ei is constant with height, consistent with the observations. However, the r ei are greater 
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than the observed (60 (im versus 25 (am). This is due partly to simulated values of n ; - that are 
smaller than the observed (Fig. 6c). 

Underestimation of n t was noticed in almost all models that participated in the ARM M- 
PACE model intercomparison (e.g. Fridlind et al. 2007). Reasons for this are not clear yet. The 
major ice forming mechanism in the CONTROL is contact-freezing of droplets. The 12 h 
averaged contact-freezing rate increases with height within the cloud layer from almost zero value 

to a value of 7 x 10 6 g kg" 1 hr" 1 at the cloud top. The formation of ice by deposition, 
condensation-freezing and immersion freezing occurs near the cloud top with an averaged rate of 

2x 10 6 g kg" 1 hr' 1 . The best quantified mechanism for ice enhancement is probably the 
shedding of ice splinters during riming, i.e., the H-M mechanism (Hallett and Mossop 1974). 
However, ice splinter production through the H-M mechanism is not significant in the simulation 

because the cloud temperature ranges from -15 °C (cloud top) to -10 °C (cloud base), colder than 

the temperature necessary for the H-M mechanism to operate (-3 °C to -8 °C). It is likely that 
other mechanisms for high ice particle concentration may be missing in the two-moment 
microphysics scheme. For example, Rangno and Hobbs (2001) argued that the fragmentation of 
delicate crystals (such as dendrites and aggregates) during crystal-crystal collisions and crystal- 
droplet collision and the shattering of some drops during freezing in free fall may play a role in 
the production of relatively high ice particle concentrations in Arctic clouds. Fridlind et al. (2007) 
claimed that two other mechanisms, formation of ice nuclei from drop evaporation residuals and 
drop freezing during evaporation, could be strong enough to account for the M-PACE 
observations. 

6. Results from sensitivity experiments 

We have shown that the CONTROL simulation reproduced most of the aircraft-observed 
cloud properties except for its underestimation of ice crystal number concentration and 
overestimation of ice crystal effective size. These two quantities have the largest uncertainties in 
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the observations (McFarquhar et al. 2007). The OneM experiment underestimated the observed 
LWC and produced a bimodal vertical structure of IWC that was not observed. Its simulated 
mixed-phase stratus glaciated earlier than in the CONTROL simulation. Four additional 
sensitivity experiments (Table 1), as described in Section 4b, are presented in this section to 
further explore the impact of microphysical processes and parameters on the simulated MPS 
clouds. Another set of sensitivity tests are used to explore the impacts of uncertainties in the 
surface fluxes, which are compared to the differences between the CONTROL and OneM 
simulations. 

6.1 Vertically integrated hydrometeor amount 

Figure 7 shows the time variability of the vertically integrated amount of each 
hydrometeor species, i.e. LWP, rain water path (RWP), ice water path (IWP), and snow water path 
(SWP) for the CONTROL, OneM, piO, NOs, SAT and the IN20 simulation. The model results are 
averaged over the entire horizontal domain in space and 30 min in time. Time-averaged values 
and standard deviations of LWP, IWP, SWP and RWP between 4 h and 12 h from all simulations 
are given in Table 2. 

A comparison among these simulations reveals the following major findings. First, a 
persistent MPS cloud layer is produced by the CONTROL, piO, and NOS simulations, which 
reaches a steady state after 3 hr, although their steady state LWP values differ (176.5±2.8, 

188.6±4.4, 142.6±6.9 g m" 2 , respectively). Second, both the OneM and SAT experiments produce 

smaller LWPs (54.2±1 5.0 and 97.9±1 2.1 g m" 2 , respectively) than the other simulations shown in 
Fig. 7a. The temporal evolutions of IWP and SWP are similar between OneM and SAT, e.g. a 
large increase near the end of simulations (Figs. 7b and c). This behavior differs markedly from 
the other simulations. Third, the IN20 experiment produces smallest time-averaged LWP 

(8.0±12.0 g m' 2 ) but largest IWP (20.9±13.2) and SWP (43.4±5.6) among the simulations, with 
LWP decreasing monotonically with time until complete dissipation at 7 h. Detailed discussions 
of these findings are given below. 
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a. Effects of spectral shape parameter of cloud ice 

The LWP derived from the MWR measurements averaged over the 12 h simulation period 
is 210 g m" 2 at Barrow, which has an approximate 20-30 g m 2 uncertainty (D. Turner 2007; 
personal communication). The steady-state LWP from the CONTROL, 176.5 g m' 2 , is 84% of the 
retrievals. The piO experiment generated more LWP (188.6 vs. 176.5 g m' 2 ) and less IWP (1.5 vs. 

4.4 g m" 2 ) and SWP (1.5 vs. 5.2 kg m" 2 ) compared to the CONTROL simulation. Decreasing the 
spectral shape parameter (p;) from 5 in the CONTROL simulation to zero in the piO experiment 
increases the phase relaxation time associated with cloud ice, i.e. making cloud ice deposition 
occur more slowly. The slower deposition to cloud ice at the expense of cloud water causes less 
IWP and SWP and more LWP in the piO experiment. Nevertheless, the difference in LWP 
between the two simulations is within the range of the uncertainties in the bulk observations that 
are presented. 

b. Effects of predicting snow number concentration 

Compared to the CONTROL, the NOS experiment produced less LWP (142.6 g m" 2 vs. 
176.5 g m" 2 ) and IWP (2.8 g m" 2 vs. 4.4 g m" 2 ) but more SWP (25.8 g m' 2 vs. 5.2 g m' 2 ). The joint 
PDF (%) of N 0s and height from the CONTROL (Fig. 8b) shows that the N 0s varies with height 

and is mostly one order of magnitude smaller than the specified constant value of 3 x 1 0 6 m" 4 in 
the NOS experiment (Fig. 8a). Other simulations with predicted N 0s produced comparable values 
of N 0s to those produced by the CONTROL. The L83 value of N 0s used in the NOS experiment 
was obtained from mid-latitude frontal cloud system. Therefore, it is not surprising that this N 0s is 
different from that predicted in Arctic mixed-phase clouds. Snapshots of the predicted N 0s (not 
shown) exhibit horizontally inhomogeneous distributions that vary with time. The significantly 
larger N 0s used in the NOS experiment resulted in stronger depositional growth of snow (i.e. 
enhanced Bergeron-Findeisen mechanism) and more significant accretion of cloud droplets by 
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snow, both contributing to less LWP and more SWP. The smaller IWP with larger N 0s is due to 
more water vapor deposited to snow and therefore less to cloud ice. 

Compared to the OneM experiment, the NOS experiment produced about the same SWP 

(~ 26 g m' 2 ) between 3 and 10 h, but less SWP after 10 h and higher LWP during most of the 
simulation period. The agreement in SWP between the OneM and NOS experiments for the 3-10 h 
period resulted from similar snow deposition rates and accretion rates of droplets by snow. The 

larger SWP in the OneM experiment after 10 h (up to 37 g m' 2 ) is related to its larger cloud ice 
deposition rates caused by the application of the water-ice saturation adjustment scheme of Lord 
et al. (1984) (Fig. 5b), and subsequently more cloud ice converted to snow. The smaller LWP of 

the OneM experiment versus the NOS experiment (54.2 vs. 142.6 g m' 2 ) is probably also a result 
from the water-ice saturation adjustment scheme utilized in the OneM experiment, as will be 
further explained below. 

c. Effects of water-ice saturation adjustment 

There are some similar results between the SAT and OneM simulations, i.e. an earlier 
decay of LWP and an increase of IWP and SWP after 9 h, compared to the other simulations 
shown in Figs. 7a-c. The relatively steady LWP between 3 h and 8 h has a magnitude of ~ 100 g 

m" 2 in the SAT experiment compared to ~ 177 g m" 2 for the CONTROL and ~ 65 g m" 2 for the 
OneM simulation. Other aspects of hydrometeor paths simulated by SAT, such as the significant 
variabilities in LWP and SWP, are more similar to those from the OneM than the CONTROL 
(Table 2). These results suggest that the application of Lord et al. (1984) saturation adjustment in 
M05 to determine the condensation of cloud water and deposition of cloud ice could significantly 
change the simulated MPS clouds and could result in an underestimation of LWP. Why did the 
SAT experiment produce less LWP than the CONTROL? It is partly because condensation rates 
near the cloud top are smaller in the SAT simulation than in the CONTROL (Fig. 9). It is also due 
to more accretion by snow resulting from a greater amount of snow in the SAT experiment. Why 
did the SAT experiment produce less IWP and more SWP than the CONTROL? Deposition rates 
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in the SAT experiment (Fig. 9) are greater than those in the CONTROL, especially near the cloud 
base and cloud top, suggesting that more cloud ice is produced from deposition of water vapor 
and subsequently more efficient conversion from cloud ice to snow. 

d. Effects of IFN number concentration 

The IN20 experiment produces greater ice crystal number concentration (~3 V 1 ) than the 

CONTROL and much closer to the aircraft observations (1-3 L' 1 ; Fig. 6c). This is attained by 
increasing the IFN number concentration to 20 times of the measured value in the IN20 
experiment. However, the increase in IFN number concentration (and hence increase in the ice 
crystal number concentration) transfers the solid, largely liquid stratus deck into a broken cloud 
system, consistent with previous modeling studies (Harrington et al. 1999; Jiang et al. 2000). As 

shown in Fig. 7e, the initially thick liquid cloud layer (LWP of 1 50 g m' 2 ) decays monotonically 
with time in the IN20 experiment and is completely gone after 7 h. This is not realistic since a 
persistent cloud layer was observed. The decay of the simulated liquid cloud layer results from the 
significantly enhanced consumption of cloud water through the Bergeron processes when the ice 
crystal number concentration is increased in the simulation. 

6.2 Precipitation and radiative flux at the surface 

Arctic clouds are linked to the hydrological cycle and oceanic processes through 
precipitation that affects fresh water input into the Arctic ocean. Surface precipitation rate is 
therefore an important parameter to be reproduced by models. Frequent, light snow events were 
reported in the ARM ground measurements during the 9-14 October period. However, as 
mentioned in Section 3b, the ARM-observed surface precipitation rate could be overestimated 
because of blowing snow (Xie et al., 2006). The temporally averaged surface precipitation rate 

from the ARM observation at Barrow was 1.7 mm day" 1 for the simulation period. Accumulated 
(horizontally-averaged) surface precipitation rates in the CRM simulations are shown in Table 3. 
The column denoted as “snow” represents the liquid water equivalent of the snow rate. The 12-h 
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averaged surface precipitation ranges from ~ 0.3 mm day" 1 to 1.2 mm day' 1 in the simulations, 
less than the ARM observations. The smallest 12-h averaged surface precipitation occurs in the 
piO and NOS experiments. However, the surface precipitation is mainly rain in the piO experiment 
and snow in the NOS experiment, respectively. The SAT and OneM experiments produce more 

surface precipitation in form of snow (~ 0.8 mm day" 1 ) than the other simulations (< 0.3 mm day" 

') except for IN20 (1.1 mm day" 1 ), consistent with their more significant ice deposition rates. 
Note that rain does not occur in the OneM experiment because the threshold for activating the 

autoconversion parameterization from cloud water to rain is 0.5 g kg' 1 , which is never reached. 
The IN20 experiment produces the largest surface precipitation among the simulations, due to the 
strong depositional growth of cloud ice and snow. 

Differences in the representation of microphysical processes affect surface radiation 
budgets through their influences on the simulated cloud microphysical and optical properties. This 
is illustrated by the temporal variations of the half-hourly and horizontally averaged downwelling 
infrared (IR) and shortwave (SW) radiative fluxes (Fig. 10). A striking feature appearing in Fig. 

10 is a dramatic reduction of the downwelling ER flux (~ 10 W m hr ) when the liquid water 
begins to decrease, i.e. after 6 h in the IN20 simulation, 10 h in the OneM simulation, and 1 1 h in 
the SAT experiment. The downwelling IR radiative fluxes in the CONTROL, piO, and NOS 

simulations differ from one another by only ~ 2 W m" 2 because these simulations all produced 
clouds that emitted as near blackbodies. The downward SW flux from the CONTROL has a 

maximal of ~ 41 W m" 2 at 5-6 h, while those from the OneM and SAT simulations differ from the 
CONTROL by up to 25 W m" 2 and 13 W m" 2 , respectively. 

6.3 Sensitivity to surface heat flux 

Figure 1 1 shows the temporal variations of LWP, IWP, SWP, and RWP when surface heat 
fluxes are changed by ±10% with either the M05 or the L83 microphysics scheme. Comparisons 
are also made against those of the CONTROL and the OneM simulations. The 4 h to 12 h 
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averages and standard deviations of the vertically integrated water amounts are given in Table 2. 
With the M05 scheme, as expected, increasing the latent heat flux and decreasing sensible heat 
flux causes generally higher hydrometeor amounts, and vice versa. However, the effects of 
changing the surface fluxes are more complicated with the one-moment scheme with larger 
variabilities, which may be related to the oscillation in the OneM simulated TKE field (Fig. 5d). 
Most importantly, it is obvious that the differences caused by the two distinct microphysics 
schemes are more significant than those due to changes in the surface heat fluxes by 1 0%. For 
example, the time-averaged LWPs from the tests utilizing the M05 scheme range between 158.2 

and 192.0 g m' 2 while those with the L83 scheme range between 51 .2 and 58.5 g m' 2 (Table 2). 

7. Summary and conclusions 

There have been few CRM simulations of boundary layer MPS clouds even though these 
clouds occur frequently in the Arctic and may potentially impact global climate and regional 
weather. In this study, a CRM has been used to simulate single-layer MPS clouds observed at the 
NS A sites during the ARM M-PACE. This CRM was implemented with both a commonly used 
one-moment microphysics scheme (Lin et al. 1983) and an advanced two-moment microphysics 
scheme (Morrison et al. 2005) as well as a state-of-the-art radiative transfer scheme (Fu and Liou 
1993). A set of simulations with different treatments of microphysical processes and different 
specifications of microphysical parameters were performed to examine the sensitivity of the 
CRM-simulated MPS clouds to cloud microphysics parameterizations. Modeled cloud fields have 
been compared to the vertical profiles of the bulk microphysical properties derived from aircraft 
measurements (McFarquhar et al. 2007), retrievals of LWP obtained from ground-based 
observations, as well as surface precipitation measurements. 

The aircraft observations suggest that the LWCs and droplet effective radii increased with 
height in the cloud layer while the droplet number concentrations, and the masses, number 
concentrations and effective radii of ice crystals were relatively constant with height. Using the 
newly implemented two-moment scheme and the observed aerosol size distribution and IFN 
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number concentration, the CONTROL simulation was able to reproduce the magnitudes and 
vertical structures of cloud liquid water, cloud ice water, droplet number concentration, and 
droplet effective radius as revealed by the aircraft observations, although it underestimated the 
number concentration of the ice crystals by an order of magnitude and overestimated the effective 
radii of the ice crystals by a factor of 2-3. With the one-moment (OneM) microphysics scheme, 
the CRM produced values of LWP and ISWP that were about 30% and 4 times, respectively, of 
those produced by the CONTROL. In addition, the vertical profile of IWC exhibited a bimodal 
distribution in contrast to the constant distribution of IWC produced with the two-moment 
approach. The deficiencies in the OneM-simulated cloud fields are largely associated with the ice- 
water saturation adjustment of Lord et al. (1984), which overestimated ice deposition rate near the 
top and base of the MPS cloud layer. The deficiencies are also closely related to the constant large 
value of snow spectra intercept parameter (n os ) used in Lin et al. (1983). Moreover, changing the 
surface latent and sensible heat fluxes, respectively, by 10%, which is an estimate of the possible 
uncertainties associated with these fluxes, caused smaller differences in the simulated cloud fields 
than those caused by application of the two- and one-moment microphysics schemes. 

When the observed IFN number concentration was used, this CRM and many other 
models that participated in a model-intercomparison project (Klein et al. 2006, Fridlinde et al. 
2007) could not reproduce the observed ice concentrations, which greatly exceeded those of ice 

nuclei (a few Zf 1 versus 0. 1 6 1' 1 ). On the other hand, the MPS clouds glaciated in the model when 
the ice concentration was initially close to the observed value, which was obtained by increasing 

the observed IFN number concentration by a factor of 20 (0.16 L" 1 to 3.2 L' 1 ). The rapid 
glaciation of cloud liquid water through the enhanced Bergeron-Findeisen process at the higher 
ice number concentration in the CRM is consistent with previous modeling studies (Harrington et 
al. 1999; Jiang et al. 2000). 

Results from a sensitivity test that used the Lord et al. (1984) water-ice saturation 
adjustment in the CONTROL simulation were more similar to those from the OneM simulation 
than the CONTROL. This indicates that the CRM -simulated MPS clouds are very sensitive to the 
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representations of cloud water condensation and cloud ice deposition. One of the major 
assumptions on which the saturation adjustment of Lord et al. (1984) is based, is that the 
saturation vapor mixing ratio (q *) is a mass-weighted average of the respective saturation values 
over liquid and ice when both cloud water and cloud ice are present. Although this assumption 
was supported by the aircraft data collected during the SHEBA/FIRE-ACE campaign (Fu and 
Elollars 2004), its utilization in models is problematic since q* depends on the model-simulated q t 
and q c and the accuracy of q t and q c prediction is influenced by other aspects of the model. Use of 
this assumption resulted in a significant underestimation of LWP during the fall season in an SCM 
simulation using SHEBA data (Yuan et al. 2006), qualitatively consistent with our findings. 

It is found that the two-moment scheme predicted much smaller values of N 0s (mostly < 

0.5 x 10 6 m' 4 ) than the constant value (3 x 10 6 m" 4 ) used in the one-moment scheme of Lin et al. 
(1983). The predicted N 0s from the two-moment approach varied with time and was both 
horizontally and vertically inhomogeneous. Using the constant larger N 0s resulted in 20% less 
LWP and 5 times more SWP. Furthermore, representing the cloud ice spectra with the exponential 
size distribution rather than the gamma distribution resulted in smaller IWP (1 .5 vs. 4.4 g m ) and 

SWP (1.5 vs. 5.2 g m ) and larger LWP (188.6 vs. 176.5 g m ), due to a slower ice deposition 
process. 

Note that the modeled results may be sensitive to some other parameters (such as the 
assumed bulk densities and fallspeeds of cloud ice and snow) involved in the two-moment 
microphysics scheme. The mechanisms for the formation of ice concentrations that greatly exceed 
those of ice nuclei in the MPS Arctic clouds and their representation in the models should be 
studied further. Additional observations and simulations are needed to further narrow down the 
uncertainties associated with these microphysical parameters because the conclusions drawn from 
this study have been based on only one case study. 
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Table 1: Description of the CRM simulations. 


CONTROL 


Microphysics 

M05 


OneM 

SAT 

IN20 


piO 

NOs 

CTR.LH+ 

CTR.LH- 

CTR.SH+ 

CTR.SH- 

1M.LH+ 

1M.LH- 

1M.SH+ 

1M.SH- 


Lin83 combined with Lord et al. (1984) and Krueger et al. (1995) 

M05 combined with Lord et al. (1984) 

same as CONTROL except IFN number concentration increased to 3.2 L' 1 
same as CONTROL except an exponential distribution is assumed for cloud ice 
same as CONTROL except for n s is diagnosed 
same as CONTROL except for increasing surface latent heat flux by 1 0% 
same as CONTROL except for decreasing surface latent heat flux by 1 0% 
same as CONTROL except for increasing surface sensible heat flux by 1 0% 
same as CONTROL except for decreasing surface sensible heat flux by 1 0% 
same as OneM except for increasing surface latent heat flux by 10% 
same as OneM except for decreasing surface latent heat flux by 10% 
same as OneM except for increasing surface sensible heat flux by 1 0% 
same as OneM except for decreasing surface sensible heat flux by 1 0% 
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Table 2. The simulated LWP, IWP, SWP, and RWP (g m" 2 ) averaged between 4 h and 12 h. The 
numbers before and after “±” are the means and standard deviations, respectively. 



LWP 

IWP 

SWP 

RWP 

CONTROL 

176.5±2.8 

4. 4+0.1 

5.2+ 0.9 

8.5+ 0.5 

OneM 

54.2±15.0 

4.7+1. 6 

28.7+4.6 

0.0+0.0 

SAT 

97.9+12.1 

4.6+0.2 

16.0+5.6 

0.3+0.0 

IN20 

8.0±12.0 

20.9+13.2 

43.4+5.6 

0.0+0.0 

piO 

188.6±4.4 

1. 5+0.4 

1. 5+0.3 

11.0+0.8 

NOs 

142.6±6.9 

2. 8+0.1 

25.8+0.7 

1.2+0. 2 

CTR.LH+ 

192.0+4.1 

4. 5 ±0.2 

6.2+1 .2 

11.9+0.8 

CTR.LH- 

158.2+4.1 

4. 4+0. 2 

4. 3 ±0.7 

5. 4+0. 7 

CTR.SH+ 

163.4+2.5 

4. 0+0.1 

6. 3 ±0.9 

6. 3 ±0.5 

CTR.SH- 

184.9+3.6 

4. 8 ±0.2 

5. 5+0.7 

9.9+0.5 

1M.LH+ 

58.0+17.1 

4.8+1 .4 

32.2+4.5 

0.0+0.0 

1M.LH- 

51.2+14.3 

3.7+1 .3 

24.6+2.4 

0.0+0.0 

1M.SH+ 

58.5+13.2 

5.2+1. 1 

27.7+1.4 

0.0+0.0 

1M.SH- 

55.2+10.3 

4.3+1. 2 

27.9+3.9 

0.0+0.0 
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Table 3. Surface precipitation rates averaged over the entire 12-h and 3-h to 12-h simulation 

periods, respectively. 


(mm/day) 

rain 

Oh-12 h 

3h-12h 

snow 

0h-12h 3h-12h 

rain plus 
0h-12h 

snow 

3h-12h 

CONTROL 

0.23 

0.18 

0.17 

0.20 

0.39 

0.38 

piO 

0.28 

0.24 

0.02 

0.02 

0.30 

0.27 

NOs 

0.03 

0.00 

0.27 

0.30 

0.30 

0.30 

SAT 

0.06 

0.00 

0.83 

0.76 

0.89 

0.76 

OneM 

0.00 

0.00 

0.76 

0.71 

0.76 

0.71 

IN20 

0.05 

0.00 

1.13 

1.28 

1.19 

1.28 
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Figure Caption 


Figure 1. (a) The r|-n c relationships represented by Eq. (5) in the text: short dashed line 
represents RL03 with a being 0.003; dot-dashed lines represent RL03± with a being 0.001 and 
0.008, respectively. Also shown are the formulations from Morrison and Grabowski (2007) (solid 
line) and Grabowski (1998) (long dashed line), (b) The corresponding p-n c relationships. See text 
for further explanation. 

Figure 2. Composite visible satellite image from the NASA Terra satellite for October 9, 
2004. The dots indicate the locations of the ARM sites at the North Slope of Alaska: Barrow, 
Oliktok Point, and Atqasuk. 

Figure 3. The upper panels show profiles of potential temperature (a), water vapor mixing 
ratio (q v ) and cloud water mixing ratio (q c ) (b) at the initial time of the simulations. The lower 
panels show profiles of the large-scale vertical velocity (c), and horizontal advective tendencies of 
temperature (d) and moisture (e), respectively. 

Figure 4. Vertical profiles of liquid water content (a) and total ice water content (b) from 
the aircraft observations (black solid lines representing the means and shadows representing plus 
and minus one standard deviation), the CONTROL simulation (red lines) and the OneM 
simulation (blue lines). Vertical profiles of ice water content (c) and snow water content (d) from 
the CONTROL simulation (red lines) and the OneM simulation (blue lines). Three lines are 
shown for each of the simulation in each panel: long dashed line 3.25 hr, short dashed line 10.25 
hr, and dot-dashed line 1 1 .75 hr. 

Figure 5. Time-height distribution of ice deposition rate (g kg' 1 hr" 1 ) sampled at 5-min 
interval from the CONTROL (a) and OneM (b) simulations, respectively. Panels (c) and (d) are 

the same as (a) and (b) except for turbulent kinetic energy (m" 2 s" 2 ) 

Figure 6. Vertical profiles of droplet number concentration (a), droplet effective radius (b), 
ice crystal number concentration (c), and ice crystal effective radius (d) from the CONTROL 
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simulation (dashed lines) and the aircraft observations (solid lines representing the means and 
shadows representing plus and minus one standard deviation). 

Figure 7. Time series of LWP (a), IWP (b), SWP (c), and RWP (d) produced by CRM 
simulations: CONTROL (solid), NOS (dots-dashed), piO (dotted), SAT (dot-dashed), and OneM 
(long-dashed with diamonds). Panel (e) represents time series of LWP (solid line), IWP (long 
dashed line), SWP (short dashed line), IWC plus SWP (dot-dashed line) and RWP (dots-dashed 
line) produce by the IN20 experiment. 

Figure 8. (a) Frequency distribution of N 0s predicted by the CONTROL simulation, (b) 
Joint PDF (%) of N 0s and height predicted by the CONTROL simulation. The contours from light 
to dark represent 0.1%, 0.5%, 1.0%, and 2.0%. 

Figure 9. Vertical profiles of cloud water condensation (solid lines) and cloud ice 
deposition (dashed lines) averaged over the 12-hr period of the CONTROL (left) and SAT (right) 
simulations, respectively. The dotted lines represent the cloud boundaries. 

Figure 10. Time-series of the half-hourly and horizontally averaged downwelling infrared 
(a) and shortwave (b) radiative flux at the surface in the CRM simulations. Panels (c) and (d) 
represent the differences between the sensitivity simulations and the CONTROL. 

Figure 1 1 . Time-series of horizontally averaged LWP (a, e), IWP (b, f), SWP (c, d), and 
RWP (d, h) from the experiments with the surface latent heat flux increased (long dashed lines) or 
decreased (short dashed lines) by 10%, or with the surface sensible heat flux increased (dot- 
dashed lines) or decreased (dotted lines) by 10%, respectively. Left panels: with the M05 scheme. 
Right panels: with the L83 scheme. The solid lines represent results from the CONTROL (left 
panels) and OneM (right panels) simulations. 
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Figure 1. (a) The r|-n c relationships represented by Eq. (5) in the text: short dashed line 
represents RL03 with a being 0.003; dot-dashed lines represent RL03± with a being 0.001 and 
0.008, respectively. Also shown are the formulations from Morrison and Grabowski (2007) (solid 
line) and Grabowski (1998) (long dashed line), (b) The corresponding p-n c relationships. See text 
for further explanation. 





Figure 2. Composite visible satellite image from the NASA Terra satellite for October 9, 2004. 
The dots indicate the locations of the ARM sites at the North Slope of Alaska: Barrow, Oliktok 
Point, and Atqasuk. 
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Figure 3. The upper panels show profiles of potential temperature (a), water vapor mixing ratio 
(q v ) and cloud water mixing ratio (q c ) (b) at the initial time of the simulations. The lower panels 
show profiles of the large-scale vertical velocity (c), and horizontal advective tendencies of tem- 
perature (d) and moisture (e), respectively. 
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Figure 4. Vertical profiles of liquid water content (a) and total ice water content (b) from the air- 
craft observations (black solid lines representing the means and shadows representing plus and 
minus one standard deviation), the CONTROL simulation (red lines) and the OneM simulation 
(blue lines). Vertical profiles of ice water content (c) and snow water content (d) from the CON- 
TROL simulation (red lines) and the OneM simulation (blue lines). Three lines are shown for each 
of the simulation in each panel: long dashed line 3.25 hr, short dashed line 10.25 hr, and dot- 
dashed line 11.75 hr. 
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Figure 5. Time-height distribution of ice deposition rate (g kg' 1 hr" 1 ) sampled at 5-min interval 
from the CONTROL (a) and OneM (b) simulations, respectively. Panels (c) and (d) are the same 

as (a) and (b) except for sub-grid turbulent kinetic energy (m" 2 s" 2 ). 
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Figure 6. Vertical profiles of droplet number concentration (a), droplet effective radius (b), ice 
crystal number concentration (c), and ice crystal effective radius (d) from the CONTROL simula- 
tion (dashed lines) and the aircraft observations (solid lines representing the means and the shad- 
ows representing plus and minus one standard deviation). 







Figure 7. Time series of LWP (a), IWP (b), SWP (c), andRWP (d) produced by CRM simulations: 
CONTROL (solid), NOS (dots-dashed), piO (dotted), SAT (dot-dashed), and OneM (long-dashed 
with diamonds). Panel (e) represents time series of LWP (solid line), IWP (long dashed line), 
SWP (short dashed line), IWC plus SWP (dot-dashed line) and RWP (dots-dashed line) produce 
by the IN20 experiment. 
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Figure 8. (a) Frequency distribution of N 0s predicted by the CONTROL simulation, (b) Joint PDF 
(%) of N 0s and height predicted by the CONTROL simulation. The contours from light to dark 
represent 0.1%, 0.5%, 1.0%, and 2.0%. 





Height (km) 



- 0.4 - 0.2 0.0 0.2 - 0.2 0.0 0.2 0.4 

Cloud Woter Condensation ond Cloud Ice Deposition (g kg -1 hr -1 ) 


Figure 9. Vertical profiles of cloud water condensation (solid lines) and cloud ice deposition 
(dashed lines) averaged over the 12 -hr period of the CONTROL (left) and SAT (right) simula- 
tions, respectively. 
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Figure 10. Time-series of half-hourly and horizontally averaged downwelling infrared (a) and 
shortwave (b) radiative fluxes at the surface in the CRM simulations. Panels (c) and (d) represent 
the differences in the infrared and shortwave fluxes, respectively, between the simulations and the 
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Figure 1 1 . Time-series of half-hourly and horizontally averaged LWP (a, e), IWP (b, f), SWP (c, 
d), and RWP (d, h) from the experiments with the surface latent heat flux increased (long dashed 
lines) or decreased (short dashed lines) by 10%, or with the surface sensible heat flux increased 
(dot-dashed lines) or decreased (dotted lines) by 10%, respectively. Left panels: with the M05 
scheme. Right panels: with the L83 scheme. The solid lines represent results from the CONTROL 
(left panels) and OneM (right panels) simulations. 





